A proposed new route to d° magnetism in semiconductors 
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Here we propose to induce magnetism in semiconductor utilizing the unique properties of the 
interstitial defect to act as the magnetic impurity within the a-PbO crystal structure. The Pbi 
interstitial generates the p-localized state with two on-site electrons to obey the Hund's rule for 
their ordering. It is demonstrated that instead of Pb interstitial other non-magnetic impurities of 
s^p^ outer shell configuration can be applied to induce d° magnetism with possibility to tune the 
local magnetic moments /is by varying a number of electrons 1 < a:: < 3. The magnetic coupling 
between such defects is found to be driven by the long-range order interactions that in combination 
with high defect solubility promises the magnetic percolation to remain above the room temperature. 
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Spintronics has recently emerged as a widely success- 
ful technology that exploits the principles of magnetism, 
but the classical metal ferromagnets are not applicable 
there. That has led to the creation of a new branch of re- 
search directed at semiconductors exhibiting magnetism 
at room temperature [l|. The first success in this regard 
was with the so-called "diluted magnetic semiconduc- 
tors" created by doping semiconductors with magnetic 
ions whose inner 3d or 4/ shells being partially filled that 
allows the ferromagnetic spin alignment [2t5| . A keen in- 
terest on ferromagnetic semiconductors also arose when 
'non-magnetic' materials were discovered to demonstrate 
magnetism 0-Q. Origin of so-called dP magnetism was 
proposed to be due to localized sp states: vacancies cre- 
ate a network of unpaired electrons and as interaction 
between such defects can demonstrate the long-range or- 
der, their magnetic coupling occurs |3]. 

Here we propose to induce d° magnetism by form- 
ing the interstitial defect Pbi (Pbi:6s^6p^ outer shell) 
within the crystal structure of the tetragonal lead ox- 
ide a-PbO. It generates the state occupied by two p- 
unpaired electrons Phi:6p^. Ordering of the localized 
electrons Pbi:6p^ obeys the first Hund's rule thus result- 
ing in formation of the stable local magnetic moment of 
2.0 /ifi. The spin-polarization energy of such state is 
Em = Eafm — i?FA/=0.235 eV, where Eafm and EpM 
are the total energies of the anti-(AFM) and ferromag- 
netic states (FM), respectively [ly]- Therefore, other- 
wise than magnetism occurs due to the partially filled p 
shell instead of the 3d or 4/ shells, Pbj works exactly 
as the magnetic impurity. It demonstrates the high on- 
site spin stability, its p electrons are localized on impu- 
rity site showing a weak perturbation with the host (an 
interaction with the host occurs through Pbi:6s^ elec- 
trons) and in addition, the interstitial defect is almost 
non-invasive to the electronic and crystal structures of 
the host. We have established that any impurity of the 
g2pi<a;<3 Q^^gj- shell embedded as the interstitial defect 
into the a-PbO crystal structure would act as magnetic. 
The p-shell occupation, pi<^<3^ can be used to control 



the spin ordering between impurities and also to tune 
the local magnetic moments: x = 2 works for 2.0 /is, 
while a; = 1 or x = 3 for 1.0 ^b- From the technologi- 
cal point of view, the layered crystal structure of a-PbO 
promises the superior advantages in achievement of the 
magnetic percolation and more importantly, a practical 
way for its control: i) the crystal structure of PbO type 
grows in polycrystalline form whose large surface area 
offers enormous potential for doping (on surface the im- 
purities can be placed on the nearest- neighbours); ii) the 
dopant solubility and long-range order interactions can 
be manipulated with the impurity atomic radius; Hi) any 
compounds of the tetragonal PbO type can be used as a 
semiconductor matrix. 

In our study we applied the generalized gradient ap- 
proximation (GGA) with the PBE parametrization 'll| 
provided by WIEN2k package for the density functional 
calculations [iJI- The supercell approach {RKmax—7) 
with sufficiently large supercell of 108-atom size (3x3x3 
array of the primitive unit cells) for single impurity and 
of 160-atom size (5x4x2) for two interacting impurities 
have been used. The Ph:5p, 5d, 6s, 6p and 0:2s, 2p elec- 
trons have been treated as the valence electrons. For 
integration of the Brillouin-zone, the Monkhorst-Pack 
scheme using a (5x5x4) k-mesh was applied. The local- 
ization of the impurity wavefunction has been addition- 
ally examined with HF applied directly to the unpaired 
electrons that allows to preserve accuracy provided by 
DFT but in the same time to correct the unpaired elec- 
trons self-interaction |13| . 

The tetragonal lead oxide a-PbO possesses the layered 
structure leading to formation of platelets upon com- 
pound growth and each platelets is considered as a sin- 
gle crystal. The layers within such crystal are held to- 
gether by the interlayer interaction of Pb:6s^ electrons 
[14| . In terms of electronic properties, these interactions 
induce a deeping of the conduction band at M* point [1^ 
and the gap shrinks as interactions enhance. GGA tends 
to overestimate the interlayer separation in the layered 
structures but in the same time it also well known to 



underestimate the band gap size. For a-PbO system it 
results in compensation effect 17] : for the lattice param- 
eters optimized with GGA the band gap is only slightly 
underestimated as 1.8 eV against the experimental value 
of 1.9 eV 1101 (application of the experimental value of 
the lattice parameters causes the band gap to shrink by 
0.22 eV). The correct band gap size is important in in- 
vestigation of the magnetic impurities: when the band 
gap size is underestimated the hole-carrying impurity or- 
bital occurs closer to the the host conduction band in 
turn forcing a delocalization of the defect tails. In or- 
der to avoid the spurious long-range order interactions 
[5| all band structure calculations have been performed 
for the lattice parameters optimized with GGA. However, 
since an incorporation of the interstitial becomes easier 
as interlayer distance increases, the interlayer distance 
has also crucial impact on the defect formation energy. 
In order to make correct evaluation of the formation en- 
ergy, the experimental value of the lattice parameters has 
been applied for those calculations. 

The interstitial defects have never been considered in 
origin of d'' magnetism because they do not belong to a 
class of most common defects in crystalline materials as 
they induce a significant perturbation of the host lattice 
resulting in their large formation energy [l^. The lay- 
ered structure of a-PbO is different, as it allows the for- 
eign atoms to squeeze between layers inducing a minimal 
lattice deformation to the host immediate neighborhood 
(see Fig. [T] (a)). As a result, the formation energy of 
Pbj is not too high. For the Pb rich conditions/vacuum 
the neutral charge state is characterized by the formation 
energy of 1.23 eV (for methods see |17|). 

The band diagram and the density of states (DOS) 
for a-PbO containing Pb^ are presented in Fig. [1] (b) 
and (c), respectively. The Pb interstitial makes a bond 
with Pb atom from the host generating the defect state 
inside the band gap. The defect-induced spin-up and 
spin-down bands (1" and l'' in Fig. [1] (b) and (c)) both 
show dominant Phi:6px+y character but demonstrate a 
different behavior. The spin-up band is filled with two 
electrons and appears just above the midgap at 0.99 eV 
+ Ey, where Ey is a top of the valence band (VB). 
In contrast, the spin-down band (l'') is unoccupied and 
appears at the edge of the conduction band (CB). In 
order to understand an asymmetry in filling of the spin- 
up and spin-down states, bonding of Pb^ with the host 
lattice has to be considered. 

The integration of the Pb interstitial into the crystal 
lattice requires excitation of its ground state 6s^6p^ to 
6s^6p^. The Phi:6s^ and Pb^iGp^ electrons participate in 
bond formation with the Pb atom from the host which 
shares its Pb:6s^ electrons (see Fig.[T](c)). The bonding 
orbital appears inside VB, while the antibonding forms 
the GB bottom (see 2" and 2'^ bands and notations of 
'bonding' and 'antibonding' at B and G panels to DOS). 
The length of the Pb^-Pb bond is fairly short, 2.9 A, that 
is indicative of a double bond formation. An out-of-plane 
displacement of 0.54 A occurs for the Pb atom involved 
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FIG. 1: (a) The Pb interstitial in the layered crystal structure 
of a-PbO (b) Band diagram: 1" and 1 are for the defect 
induced bands, while 2" and 2"* are the antibonding orbitals 
of the Pbi-Pb bond, (c) Total and partial density of states for 
Pbi and host Pb atoms ('top' is referred to Pb atom forming 
the Pbi-Pb bond, while 'bottom' for Pb atoms in the bottom 
layer). The 'bonding' and 'antibonding' in panels B and C 
are reffered to the Pbi-Pb bond. 



in bonding. Because participation of the host Pb:6s^ 
orbitals in bonding interferes with interlayer interactions, 
it alters a behavior of the conduction bands near the GB 
top : in analogy with the band behavior in a single layer 
p^ . the band deeping at the M* point is suppressed. 

Therefore, out of three Pbi:6p'^ electrons found in the 
excited state {6s^6p^), only one Phi:6pl participates in 
bonding while the two others Phf.Gpx+y are left at the 



defect site. The Hund's rule dictates both unpaired elec- 
trons to occupy the 1" spin-up band. The empty l'' spin- 
down band is pushed up to CB that causes considerably 
large spin-exchange splitting of order 0.523 eV. The fer- 
romagnetically ordered Phi:6px+y electrons generate the 
local magnetic moment of 2.0 fis- The difference in to- 
tal energy between the AFM and FM states is found 
to be Em — Eafm — i^FM =0.235 eV thus promising 
the FM state to be stable well above the room tempera- 
ture. An application of the HF approach to the unpaired 
electrons induces further stabilization of the FM state to 
£'M=0.490 eV due to stronger on-site localization of the 
impurity wave function (the 1" orbital is shifted towards 
the valence band by 0.5 eV thus enhancing the splitting of 
the 1" and 1"^ orbitals). On the other hand, reduction of 
the interlayer distance (the experimental lattice parame- 
ters are applied) in opposite causes delocalization of the 
wave function due to stronger hybridization of impurity 
state with the host lattice. 

The 6s^6p^ electronic configuration is a key point for 
the Pb interstitial to act as a magnetic impurity: filled s^ 
outer shell is required for bonding with the host lattice 
while partially filled p shell contributes in development 
of the local magnetic moment. Following this principle, 
other chemical elements possessing s^p^-^~^ outer shell 
can be applied to induce dP magnetism. Indeed, several 
examined impurities have shown a formation of the local 
magnetic moment, impurity of x = 1 or a; == 3 induces 
the local magnetic moments of 1 /is, while x = 2 pro- 
duces 2 fiB- The stability of the spin-polarized state is 
observed to grow substantially with reduction of the im- 
purity atomic radius. 

Our studies of localization of the electron density at the 
defect site has revealed a duality in its behavior. Since 
occupied band 1" is located deep inside the band gap, 
the Phi:6p'^^y electrons shows strong localization on the 
defect site thus allowing to form the stable local magnetic 
moment of 2/iB (see Fig. [2] (a)). However, the long defect 
tails appear as well but due to bonding/hybridization of 
the Pbi:6s^ and Phi'.Gpl electrons from the Pb intersti- 
tial and Pb:6s and Pb:6p, 0:2p electrons from the host 
(see partial DOS in Fig. [T](c)). As a result, the electron 
density around the defect site is spin-polarized showing 
an anisotropy in its distribution: the defect tails are sym- 
metrically polarized relative to Pbi thus contributing in 
stabilization of the ferromagnetic state. 

Since extension of the wave function tails to the host 
lattice seems to be of the long-range order, it is expected 
to induce effective defect-defect interactions. In this re- 
spect we have investigated an interaction of Pb^ with the 
O vacancy (Vb) which formation is also favored at the 
Pb-rich/0-poor limit (oxygen deficiency). In its charge 
neutral state, Vq is occupied by two electrons showing the 
strong localization at the vacancy site [ITj . Since the vir- 
tual hopping is allowed between the states formed by Pb^ 
and Vb, they can work as compensating centers to each 
other. For the defects appearing at the nearest-neighbour 
sites, the FM ordering of the localized spins at Pb^ is 
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FIG. 2: Schematic band diagrams demonstrate a formation 
of FM (Phi) and AFM (Pbi-Vb pair) states. The energy 
values are appearance of the defect levels relative to Ev- The 
electron density map is plotted for isovalues of ±0.005 e/A'^ 
with help of Xcrysden. (a) FM state at the Pbi site. The 
density map demonstrates the spin density calculated for the 
energy range (1.0 eV-|--Bv)±0.15 eV. (b) An interaction of Pbi 
with Vb destroys the FM state at the Phi site. The electron 
density is plotted for the energy range (0.9 eV+iJv)±0.4 eV. 



found to be destroyed due to the strong interaction be- 
tween defects to be accompanied by hybridization. The 
Phf.dpx+y electrons from the interstitial and the Pb:6s 
and Ph:6p and 0.2p electrons from the host lattice, all 
have been found to contribute equally in formation of 
both defect states (see manifestation of the electron den- 
sity in the inter-defect region in Fig. H^b)). The strong 
interactions between defects contribute in lowering of the 
formation energy of defect pair by 0.84 eV in comparison 
to the non-interacting defects. With defect separation 
the electronic interactions are suppressed leading to re- 
covery of the FM state: for distance of 6.2 A between the 
defects the local magnetic moment is 1.8 fiB- However, 
because interaction of Pbi with Vq do not vanish com- 
pletely, the magnitude of the spin polarization energy is 



found to be considerably low Em = Eafm — -E_fa/=0.06 
eV. Therefore, the electronic interactions between Pb^ 
and Vb demonstrate the long-range order due to the long 
tails of the Pb^ interstitial. 

For two interacting interstitial defects, as both defects 
contribute with their tails into the long-range order inter- 
actions, the magnetic coupling between them is expected 
to sustain over the larger distance. The Gp^+y state is 
exactly half filled and, therefore, for two interacting im- 
purities the virtual hopping is allowed only in the AFM 
state making it the ground state Q . For two Pb intersti- 
tials placed on distance 4.1 A we found that -Em=-0.96 
eV while for distance 12.5 A it drops down to -0.0056 eV. 
For the C interstitial the interactions are lower as Em—- 
0.38 eV for distance 4.1 A and Sm=-0.0023 eV for 12.5 
A. An application of HF for the localized electrons sup- 
presses Em almost twice. The AFM interaction between 
impurities can be switched to FM for the impurity of p^ 
occupation of outer shell (for example. In or Ga work for 
FM state) or for p^ occupation by choosing the charged 
state (1+). 

In terms of the formation energies of the interstitials, 
it drastically drops down as the atomic radius decreases. 
Thus, if the formation energy of the Pb interstitial is 1.23 
eV for the Pb-rich/vacuum conditions, it is reduced al- 
most to zero for Ge, Si impurities and already turns to be 
negative for the C and O impurity (the formation energy 
of the O interstital is -0.26 eV) thus implying the higher 
solubility limit with prospects of spontaneous defect ap- 
pearance. Therefore, it is found that suppression in the 
magnetic coupling between impurities to occur with re- 
duction of its atomic radius is compensated by simulta- 
neous shift of the impurity solubility limit to the higher 
defect concentration. For defects appearing on surface, 
the defect concentration can reach a number of sites avail- 
able for doping {^ 10^^ cm~^) as the formation energy is 
reduced further down (for example for the Pb interstitial 
it drops down by ^ 1.0 eV such as the solubility limit is 
10^" cm~^). 

In summary, we propose to induce the p local orbital 



magnetism by doping of a-PbO semiconductor with the 
non-magnetic impurities to appear in layered structure 
of a-PbO as the interstitial defects. To create conditions 
for partially filled p shell to act as 3d or 4/ shells 
of the magnetic ions, i.e. obeying the Hund's rule 
for the on-site ordering of the unpaired electrons, an 
impurity of specific s^p^-^-^ outer shell is required. 
The combination of the outer shell with the PbO crystal 
structure is unique because it allows doping to be almost 
'non-invasive' to the host lattice and its electronic prop- 
erties. In particular, bonding between impurity and the 
host involves only their s^ outer shell electrons thereby 
preserving the original electronic configuration of the p 
shells. In this case, the partially filled p^-^-^-shell of 
impurity generates the localized spins on-site of defect. 
The main advantage of the defect-induced magnetism 
over the diluted magnetic semiconductors is duality 
in its state localization. Thus, the p-localized state 
formed on site of the Pb interstitial has the localized 
nature to form the stable local magnetic moment, but 
impurity-host hybridization results in the extended de- 
fect tails. The manifesting long-range order interactions 
between defects in combination with their high solu- 
bility create conditions for magnetic percolation to occur. 
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